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Based on the assumption that thin layers of fee CrAs can be sandwiched lattice-matched in 
between GaAs we perform a theoretical study of spin filtering in [1,0,0] GaAs/CrAs/GaAs het- 
erostructures. The electronic structure of fee GaAs, CrAs, and GaAs/CrAs heterostructures is 
obtained within LMTO local spin-density functional theory. Relevant segments of the computed 
band structure are mapped onto an effective nearest-neighbor sp 3 d 5 s* tight -binding model from 
which, together with the computed band offset, the effective Hamiltonian of the heterostructure 
is constructed. The current -voltage characteristics are computed within a non-equilibrium Green's 
function formalism and for several temperatures and thicknesses of the CrAs layer. The results from 
this investigation indicate that high spin polarization should be feasible with ideal GaAs-CrAs-GaAs 
heterostructures and persist even at room temperature. 
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I. INTRODUCTION 

The identification and design of high-efficiency all- 
semiconductor spin filtering devices which operate at 
room temperature and at zero external magneti c fie ld 
are of profound interest for spintronic applications!^ In 
spite of some success, for example regarding improved 
spin injection from ferromagnets into Si, progress has 
been moderate. Dilute magnetic semiconductors, such 
as GaMnAs, which have been grown most successfully 
with good quality of interfaces onto nonmagnetic fee 
semiconductors, generall y, h ave critical temperatures far 
below room temperature.- 3 Growing high-quality GaM- 
nAs layers within heterostructures proves to be more 
difficult than in the bulk. Most likely, due to struc- 
tural problems in the form of unwanted defects lead- 
ing to free-carrier compensation, an inevitable disor- 
der in ternary alloys, and limitations to the hole dop- 
ing levels in the semiconductor contact layers, in exper- 
iment, ferromagnetic order appears to be difficult to es- 
tablish and to maintain under bias in AlGaAs/GaMnAs 
heterostructures!^ Signatures of (weak) quantum con- 
finement effects associated with GaMnAs quantum wells 
in GaAlAs/GaMnAs single- and double barrier het- 
erostructures have been reported in the literature.^ Us- 
ing one Gao.96Mno.04As contact layer on an asymmet- 
ric GaAlAs/GaMnAs heterostructure magnetization- 
dependent negative-differential-resistance has been 
observed.^ A similar TMR experiment has been per- 
formed recently, in which resonant tunneling in non- 
magnetic AlGaAs/GaAs/AlGaAs was used to control the 
hole current P 

MAs or MSb compounds, with M denoting a tran- 
sition metal, such as V, Cr, Mn in the bulk are not 



stable in the fee phase. Even in thin layers on an fee 
substrate, such as GaAs, they appear to be difficult to 
grow in this phasePtHl Some noteworthy exceptions are 
the reports of the experimental realization of fee MnAs 
quantum dots on GaAs a nd th e growth of thin layers of 
CrAs on GaAs substrates! 13 * 14 * Recent experiments have 
suggested that CrAs is a halfmctallic ferromagnet which, 
when grown epitaxially in the fee structure on top of 
GaAs, displays ferromagnetic behavior well above room 
temperatureP^Hlfil Ab-initio studies of the system have 
led to the prediction that the C urie t emperature of fee 
CrAs may be as high as 1000 KpEHU The electronic 
structure of GaAs / CrAs heterostructures has been stud- 
ied by Bengone et alW^. They also studied transport 
properties through GaAs/CrAs/GaAs trilayers. R has 
to be noted at this point that density functional theory 
(DFT) based calculations, as well as some experimental 
resultJ^JHI indicate that the fec-structure of thin film 
CrAs is energetically unstable. R has therefore been ar- 
gued that the experimentally observed ferromagnetic be- 
havior may be caused by magnetic defects near the het- 
erointerface and that no half-metallicity may be present 
in thin layers of CrAs grown on GaAs substrates at all. 

In the recent past we have investigated the effect 
of electronic correlations upon the half-metallicity of 
stacked short period (CrAs)f/(GaAs)f (£ < 3) superlat- 
tices along [001]. Our results indicate that the minority 
spin half-metallic gap is suppressed by local correlations 
at finite temperatures and continuously shrinks on 
increasing the heterostructure period. As a consequence, 
at the Fermi level, the polarization is significantly 
reduced, while dynamic correlations produce only a 
small deviation in magnetization^. Essentially the 
same effect was found for defect-free digital magnetic 
heterostructures £-doped with Cr and MrP2. In addition 
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both studies indicate that the minority spin highest 
valence states below the Fermi level are localized more 
in the GaAs layers while the lowest conduction band 
states have a many-body origin derived from CrAs. 
Therefore independent of the presence of electronic 
correlations in these hctcrostructures holes and electrons 
may remain separated among different layers which 
may be detected in photo-absorption measurements. 
From these combined density functional and many-body 
calculations^^ briefly described above, it is not 
possible to establish whether the many-body-induced 
states within the minority spin channel are conducting 
or not. It is of great interest to study transport through 
such heterostructures. 



The aim of the present paper is to explore spin filtering 
in GaAs/ single-layer-CrAs/GaAs heterostructures un- 
der the assumption of perfect fee lattice matching. As 
will be shown below half-metallicity is not a necessary 
ingredient for efficient spin filtering. Essential are the in- 
terfacial properties and the spin-selective band matching 
(bias-dependent "band alignment") between the CrAs 
and GaAs layers. Within the current work we determine 
the band matching by aligning the bands according to 
the band offset obtained within the local spin-density ap- 
proximation (LSDA) results. The interfacial properties 
are modeled ad hoc by keeping constant the As onsite 
energies appearing in our tight-binding model through- 
out the device. We demonstrate in the current work that 
ideal GaAs /CrAs/GaAs heterostructures should func- 
tion as room temperature spin filters regardless of half- 
metallicity. 



Our theoretical approach consists of three main parts: 
In a first step we determine the spin-dependent band 
structure of lattice matched fee bulk CrAs and a [1,0,0] 
(GaAs)6 / (CrAs) 6 supercell using the LMTO method. 
This method is also used to calculate the electronic struc- 
ture of GaAs. From these calculations we determine the 
band offset between the two materials, as well as the 
spin-dependent bulk electronic structure of CrAs. We 
then map the electronic structure of the bulk materi- 
als onto an effective 20-orbital sp 3 d 5 s* nearest neighbor 
tight-binding (TB) model. In the final step we calcu- 
late the current-voltage (I-V) characteristics within the 
non-equilibrium Green's function formalism. 



This article is structured as follows. In Sec. |TT] we 
present the theoretical approach and the results for the 
electronic structure. In Sec. |III| we describe the mapping 
procedure onto the tight-binding model. The transport 
model and results for several layer thicknesses and tem- 
peratures are discussed in Sec. |IV| Summary and con- 
clusions are given in Sec. [Vj 



II. LMTO ELECTRONIC STRUCTURE 
CALCULATIONS AND THE CRAS/GAAS [1,0,0] 
BAND OFFSET 

The electronic structures of bulk zinc-blend (ZB) 
GaAs, bulk ZB CrAs as well as of a (GaAs) 6 (CrAs) 6 
supercell have been determined employing the L MTO- 
ASA code as developed by Jepsen and Andersen! 24 ! 25 ! 
This code has been employed previously to explore the 
electronic structure of bulk ZB CrAs and thin-layer fee 
GaAs/CrAs superlattices. In particular the half-metallic 
behavior as a function of s uperl attice period and lattice 
constants was investigated! 22 ! 26 ! Details of this approach 
can be found in the literature.^ 

The electronic structure model is based on LSDA omit- 
ting spin-orbit interaction and corrections for strong cor- 
relations as provided, for example, by the dynamic-mean- 
field-theory (DMFTp-3 or variational cluster approach 
(VCA}2SES1. The omission of these corrections clearly 
leaves open the question as to the accuracy of the results 
fo r fee CrAs. However, it has been shown by Chioncel et 
al! 22 ! 26 ! that the inclusion of correlations does not affect 
the magnetization. LSDA or its gradient-corrected ap- 
proximation (GGA) produce band gaps that are typically 
at least 30% smaller than the experimental values for al- 
most all semiconductors and insulators. In the case of 
half-metals one may speculate that this problem is not se- 
rious, as the dielectric response of half-metals is of metal- 
lic type. State-of-the art GW calculations show the ten- 
dency to improve on the gap problem for semiconductors, 
however, applied to half-metallic Lao.7Sro.3Mn03 they 
over-estimate the half-metallic gapP3 The well-known 
GaAs band gap of 1.52 eV at low temperature is strongly 
underestimated by the present spin-DFT method, pre- 
dicting a value of about 0.35 eV, while the overall features 
in the vicinity of the band gap are reproduced reasonably 
well. Since we do not study transport across the main 
energy gap, such as in Zener tunneling, in what follows 
we consider two distinct cases: (i) we leave the electronic 
structure unchanged, i.e. we use the GaAs band struc- 
ture with the underestimated band gap, and (ii) we scis- 
sor the GaAs electronic structure to the experimentally 
observed value of 1.52 eV. For want of better knowledge, 
the electronic structure of majority and minority spin 
fee CrAs was not scissored. Previous results showed that 
the detailed nature of states around Ep in bulk CrAs 
is changed when different lattice constants or models for 
electronic correlations are considered.^ Furthermore, we 
consider n-doped GaAs such that the spin-orbit effects, 
at least in GaAs, can be neglected. Here clearly the fi- 
nal answer as to the accuracy of the electronic structure 
model can only be given by comparison to experiment. 
With the present work, in fact, we hope to kindle interest 
in further experimental studies of this material combina- 
tion. 

Since the GaAs band structure is well known and, ex- 
cept for the band gap, is fairly well reproduced within 
the LMTO-ASA code, here we mainly discuss the results 
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obtained for bulk fee CrAs with the GaAs lattice con- 
stant ai — 5.65 A. The electronic structure for majority 
and minority carriers is shown, respectively, in Figs, [l] 
and [2] (red solid lines). For the minority carriers a gap of 
about 1.8 eV is predicted. A more detailed analysis of the 
electronic structure shows the origin of the valence band 
edge in the As-p - Cr-d hybrid orbitals while conduction 
band edge sta tes ar e dominated by the Cr-d orbitals, see 
Chioncel et alP^l. 

Let us emphasize an interesting point: We observe 
from the ab-initio calculations that the lattice constant 
chosen for bulk fee CrAs d eterm ines whether CrAs is 
half-metallic or not, see alsd 22 l 26 l Since the bulk band 
structures of GaAs and CrAs are shifted relative to each 
other by the band offset, as discussed below, we obtain 
that half-metallicity of the CrAs layer is not compulsory 
for the realization of a spin filter. A necessary ingredient 
is a high degree of spin-dependence in the CrAs band 
structure. 

In the spirit of core level spectroscopy, the band align- 
ment (band-offset) at the [1,0,0] GaAs/CrAs interface 
was determined using the average energy of the low-lying 
As-s bands of the two bulk materials for a reference state. 
We computed the center of energy of these bands for bulk 
fee CrAs Ac, bulk fee GaAs Aq (relative to the valence 
band edge), as well as for the As atoms embedded in the 
Ga and Cr environment in a (GaAs)e(CrAs)6 supercell 
Acg and Acc (relative to the Fermi energy). The center 
of energy of one particular band is obtained by calculat- 
ing the fat-band weights with the help of the ab-initio 
LMTO-ASA code and using them as a weight when per- 
forming the average over the first Brillouin zone. The 
band offset A is calculated according tcP^ 



significantly change the spectral overlap of the lowest 
GaAs conduction band with the CrAs bands near the T 
point. The situation is rather different for minority car- 
riers, however, here the overlap of the lowest GaAs con- 
duction band in both cases occurs with fairly flat CrAs 
bands of low mobility. 
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Figure 1: (Color online) Majority spin CrAs ab-initio band 
structure (red solid line), scissored, and un-scissored GaAs 
band structure (blue dashed line and green dash-dotted line) . 
The Fermi energy Ep — 0.01 eV above the conduction band 
minimum of GaAs is indicated by the horizontal lines (solid: 
scissored GaAs, dashed: un-scissored GaAs). 



A = Aq — Ac + A GG - A 



cc- 



(1) 



With this procedure we obtain Ai = 0.60 eV and 
A2 = 0.55 eV, respectively, for majority and minority 
spin. Here it should be recalled that the determination of 
the band-offset is based on a (GaAs)e(CrAs)g supercell. 
A reasonable error estimate is difficult, since it involves 
not only the procedure of finding the center of energy for 
the As-s band but also the underlying electronic struc- 
ture model. In the present calculations we have used the 
LSDA-ASA (atomic sphere approximation) method and 
the potentials were converged to better than 10 -6 Ry, 
with the same accuracy for the center of the bands and 
related quantities. 

In Figs, [l] and [2] we show, on top of the electronic 
structure of majority- and minority-spin CrAs, the elec- 
tronic structure of GaAs for case (i) (un-scissored gap) 
and case (ii) (scissored gap) taking into account the com- 
puted band offset. It can now clearly be seen that carriers 
injected from n-GaAs near the T-point do not reach the 
main energy gap region of bulk CrAs. Even in case of p- 
doping of the GaAs contact layers, this region can only be 
reached under an applied bias exceeding approximately 
0.3 V, as can be observed from an inspection of Fig. [2] 

Moreover, scissoring of the GaAs band gap does not 
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Figure 2: (Color online) Minority spin CrAs ab-initio band 
structure (red solid line) and scissored as well as un-scissored 
GaAs band structure (blue dashed line and green dash-dotted 
line). The Fermi energy Ef = 0.01 eV above the conduction 
band minimum of GaAs is indicated by the horizontal lines 
(solid: scissored GaAs, dashed: un-scissored GaAs). 
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III. EFFECTIVE SP 3 D 5 S* TIGHT-BINDING 
MODEL 

In the second step the ab-initio electronic structure 
£ n SDA °f GaAs as well as of majority- and minority-spin 
CrAs individually is mapped onto an effective sp 3 d 5 s* 
nearest neighbor TB modelP^U 

In principle, this step can be avoided if the LSDA wave 
functions were used to express the transmission function 
of the heterostructure, as for example proposed within 
the SIESTA DFT approachP The trade-off of an ap- 
proach which is based on the L(S)DA one-particle wave 
functions, however, is that its validity ad hoc is question- 
able, since the wave functions used in the Kohn-Sham 
variational principle do not allow a direct physical inter- 
pretation so that their connection to the S-matrix is not 
obvious. The use of ground state wave-functions would 
definitely limit one to the linear-response regime, since 
the transmission function would be obtained for zero ex- 
ternal bias only. In general LDA bulk band-structure 
calculations are flawed regards to producing the correct 
energy gap. While for bulk GaAs it is known which scis- 
soring procedure is required to fix this problem, for fee 
CrAs it is not. For an ab-initio supercell calculation 
for GaAs/CraAs heterostructures a simple "scissor op- 
eration" does not work at all. These difficulties have 
convinced us to follow a mapping (" downfolding" ) ap- 
proach from the ab-initio band structure calculation to 
an empirical tight-binding model. 

The benefits from downfolding onto an tight -binding 
model are that, for the energy window of interest, an 
effective single-particle Hamiltonian for the individual 
bulk systems, as well as the bulk heterostructure in their 
respective ground state is constructed. For each bulk 
system an individual scissor operation can be employed 
if deemed necessary. Subsequently, the Hamiltonian of a 
given heterostructure can be put together layer by layer. 
Note that the calculation procedure for the band align- 
ment (here between GaAs and CrAs) is independent of 
scissoring. Moreover, an applied electric bias and bias- 
dependent mean-field corrections can be added in self- 
consistent fashi on, a s utilized by some of us recently for 
similar systems! 36 * 37 * This convenient layer-by-layer con- 
struction can be employed in the construction of the 
non-equilibrium Green's function components also and 
has been used in the calculations below. A further ad- 
vantage of the approach chosen within this work is that 
we know the /c-dependence of the bulk TB Hamiltoni- 
ans analytically via the structure factors. Hence, the 
influence of small deviations in the TB binding param- 
eters can be investigated in a systematic manner. On 
the other hand, the tight-binding parameters cannot be 
determined uniquely and must be understood as fitting 
parameters only. Moreover, there is no clear prescrip- 
tion of how to interpolate them at a heterointerface. A 
detailed study of the benefits and caveats of such an ap- 
proach will be given in the near future.^ 

The formal process of downfolding is executed using a 



genetic algorithm as implemented in Matlab and mini- 
mizing the cost functional 

K (0 = /$>n(fc) [ £ LSDA {k) _ £ TB M ^ . (2 ) 
y nk 

Here, a n (k) are normalized weights where n is the band 
index and k the wave vector, £ denotes the set of 31 in- 
dependent TB parameters and ej s (fc,£) is the TB band 
structure as a function of k and £. The weights a n (k) 
are used to restrict and/or focus the fit to the part of 
the band structure which contributes to charge transport. 
This ensures that the energy bands are well represented 
by the TB fit and no " spurious bands" appear inside this 
energy window. 

The fit procedure is executed in two steps: In a first 
step we fit the majority and minority spin band struc- 
ture of CrAs independently without any further restric- 
tions. In a second step we fit the GaAs band structure 
under the constraint that all As onsite energies have to 
have the same value as in CrAs. Let us briefly discuss 
the implications of such a fitting procedure: First of all, 
we obtain two different sets of TB parameters for GaAs 
since we fit the two CrAs band structures independently 
and then restrict the GaAs parameters depending on spin 
orientation. Had we included the spin-orbit interaction 
into our model, two distinct sets of TB parameters for 
GaAs would come more natural, however, in any case 
the values obtained must be considered best fits under 
given constraints. Such a constraint on the GaAs TB 
parameters is motivated as follows: the As onsite ener- 
gies have to be equal in both materials in order to avoid 
the question of how to choose the onsite energies of As 
near the GaAs/CrAs interface, in particular, the inter- 
face As sandwiched between a Cr and a Ga cation. Since 
such a fitting procedure in general will not yield unique 
TB parameters, this constraint also helps to limit the 
number of possible parameter sets. 

All in all we perform the downfolding process for four 
different combinations: (A) majority-spin CrAs and scis- 
sored GaAs, (B) majority-spin CrAs and un-scissored 
GaAs, (C) minority-spin CrAs and scissored GaAs, and 
(D) minority-spin CrAs and un-scissored GaAs. The 
TB parameters which were identified as optimal for each 
of the four cases are listed in the Appendix. The fits 
for GaAs together with the ab-initio band structure and 
the fits of majority- and minority-spin CrAs are given 
in Figs. [3J |4j [5] and [6] In view of the fact that the com- 
puted ab-initio electronic structure can at best capture 
the overall features of the actual electronic band struc- 
ture the obtained TB fits are highly satisfactory. 

In principle, the spin-orbit interaction can be included 
in the TB model following the work of ChadiP^ However, 
since the LMTO-ASA code itself currently does not fea- 
ture spin-orbit interactions, its implementation at the 
TB level would require the introduction of further (and 
somewhat arbitrary) parameters into our model, in par- 
ticular for CrAs. Furthermore, we shall focus on trans- 
port with n-doped GaAs buffer layers so that a detailed 
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account of the spin-orbit interaction in the GaAs elec- 
tronic structure will not really be important here. 
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Figure 3: (Color online) Scenario (A): majority spin CrAs 
ab-initio band structure (red solid line), CrAs TB-fit (blue 
dashed line) , scissored GaAs band structure (green solid line) 
and GaAs TB fit (black dashed line). The Fermi energy 
Ef = 0.01 eV above the conduction band minimum of GaAs 
is indicated by the horizontal solid line. 
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Figure 4: (Color online) Scenario (B): majority spin CrAs 
ab-initio band structure (red solid line), CrAs TB-fit (blue 
dashed line), un-scissored GaAs band structure (green solid 
line) and GaAs TB fit (black dashed line). The Fermi energy 
Ef = 0.01 eV above the conduction band minimum of GaAs 
is indicated by the horizontal solid line. 
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Figure 5: (Color online) Scenario (C): minority spin CrAs 
ab-initio band structure (red solid line), CrAs TB-fit (blue 
dashed line), scissored GaAs band structure (green solid line) 
and GaAs TB fit (black dashed line). The Fermi energy 
Ef = 0.01 eV above the conduction band minimum of GaAs 
is indicated by the horizontal solid line. 
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IV. STEADY-STATE TRANSPORT 

The I-V-characteristics of the heterostructure is cal- 
culated within a non-equilibrium Green's function ap- 
proach which has been adapted from our recent study of 
GaMnAs-based heterostruct ures, to w hich we refer for where Lbz 
further details and references! 36 * 37 * 41 ^ ^ 



Figure 6: (Color online) Scenario (D): minority spin CrAs 
ab-initio band structure (red solid line), CrAs TB-fit (blue 
dashed line), un-scissored GaAs band structure (green solid 
line) and GaAs TB fit (black dashed line). The Fermi energy 
Ep = 0.01 eV above the conduction band minimum of GaAs 
is indicated by the horizontal solid line. 



The Hamiltonian of the GaAs/CrAs/GaAs het- 
erostructure is obtained by performing a partial Wannier 
transformation from the wave vector k to (x, fey), where 
x denotes the [1,0,0] growth direction of the crystal and 
fen is the in-plane (parallel) k-vectoi^S 
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— and 7 labels the layer. 
The resulting one-particle Hamiltonian is of the single- 
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particle form 

2,(7(7' 

+ E < -'( fc iiH + i, CT ( fc iiK-'( fc ii) + h - c - I ( 4 ) 

i,<7<r' 

with c[ a {k\\) ( Ci iCT (fc||)) denoting the creation (annihi- 
lation) operator for site i and orbital a. (fey) and 

tfl,{k\\) are onsite and hopping matrix elements. The 
semi-infinite GaAs leads are taken into account by eval- 
uating the associated self-energies and feeding them into 
the system's Dyson equation!^ The surface Green's func- 
tions are obtained with the help of an algorithm sug- 
gested by Sancho et alP^. The transmission function 
T(E, fcy) for total energy E and in-plane momentum fcy 
is calculated via 



T{E,k\\) = Tr [T R G R T L G A ] 



(5) 



Here, G R I A are the system's retarded (R) and advanced 
(A) Green's functions, T^m are the coupling functions 
to the left (L) and right (R) GaAs leads and Tr [•] is 
the trace operation. We then compute the steady-state 
current j(V a ) assuming local thermal equilibrium among 
the electrons injected from a particular contact using the 
standard expression from stationary scattering theory^ 



j(V a ) 



2e 



£ / dET(E,k {l )[f L (E)-f R (E)], 



(6) 



with e, h and Jl/r denoting, respectively, the elemen- 
tary charge, Planck's constant, and the Fermi-Dirac dis- 
tribution function for the left and right electric contact. 
V a enters Eq. (|6| in two places: the transmission func- 
tion T(E, fcy) and the difference in the quasi-Fermi levels 
between left and right contact. In order to cut com- 
putational cost, we assume a linear voltage drop from 
the left to the right lead across the simulated structure. 
This implies a somewhat artificial relationship between 
the electr ic field across the structure and the applied 
bias! 44 * 45 J In principle, both an effective single-particle 
potential, an effective exchange splitting, and a self- 
consistent treatment of charge injection can be imple- 
mented into the present model.^ However, while provid- 
ing a significant reduction in computation time, omis- 
sion of self-consistency does not significantly reduce the 
quality of our results in view of other approximations 
made. Note that the Fermi-Dirac distributions of the 
GaAs contacts provide the sole temperature dependence 
in the current model since a temperature dependence of 
the electronic structure is not considered here. 

In what follows we present results for the I-V char- 
acteristics of GaAs/(CrAs)f/GaAs heterostructures for 
I = 4,6,8,10. While thin layers of fee CrAs may be 
easier to realize experimentally, thicker layers thereof are 
more realistically described within our approach. The 



free carrier density in the n-doped GaAs regions is about 
4.5 x 10 17 cm- 3 at T = 300 K (7.9 x 10 16 cm" 3 at T = 
K), with the quasi-Fermi level held constant at 10 meV 
above the conduction band edge. The applied bias was 
varied between zero and 0.2 V. Results for scissored and 
un-scissored GaAs, respectively, and I = 10 are shown 
in Figs. [7] and [8] (mind the semi-logarithmic plot). The 
overall features of the I-V characteristics agree for both 
cases: the majority current density clearly dominates the 
minority current density and, in most bias regions, by 
several orders of magnitude. However, this effect is more 
pronounced for the scissored GaAs model. While the ma- 
jority current is rather insensitive to scissoring, the mi- 
nority current density is not (we believe that the small 
oscillations for the 77 K minority case near 1.5 V in Fig. 
[7] are of numerical origin) . The reason is found when in- 
specting Figs, [l] and [2] It shows that, for un-scissored 
GaAs at low applied bias, there is a resonance between 
the GaAs conduction band minimum (dashed line) and 
CrAs associated bands near the T-point. Near the F- 
point these bands are rather flat and so the group veloc- 
ity is almost zero. Under moderate bias, however, these 
bands are moved further into resonance (to regions with 
higher group velocity) with the conduction band of GaAs 
at the emitter side and the minority current initially rises 
steeply. 




Figure 7: (Color online) Spin-resolved current-voltage charac- 
teristics for GaAs/ (CrAs) 10/GaAs and case (A) and (C) (scis- 
sored GaAs) for different temperatures (T = 4.2 K, T = 77 K, 
and T = 300 K). The majority spin current clearly dominates 
the minority spin current for all voltages and temperatures. 

The current spin polarization P(V a ) as a function of 
applied voltage V a is defined as 



p(y a ) 



Jmin (Va) 



Jmaj (Va) + Jmin (Va) 



(7) 



Here, j ma j/ m in refers to the majority and minority spin 
current density, respectively. In Figs. [9] and |10| respec- 
tively, we display the computed current spin polarization 
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Figure 8: (Color online) Spin-resolved current- volt age char- 
acteristics for GaAs/(CrAs)io/GaAs and case (B) and (D) 
(un-scissored GaAs) for different temperatures (T = 4.2 K, 
T = 77 K, and T = 300 K). The majority spin current clearly 
dominates the minority spin current for all voltages and tem- 
peratures. 



for scissored and un-scissored GaAs and the three differ- 
ent temperatures discussed above. 
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Figure 9: (Color online) Current spin polarization P(V a ) for 
scissored GaAs for different temperatures (T = 4.2 K, T = 
77 K, and T = 300 K). 



Figure 10: (Color online) Current spin polarization P(V a ) for 
un-scissored GaAs for different temperatures (T = 4.2 K, 
T = 77K, andT = 300 K). 



slightly changed if our simulations are performed under 
equal-electric-field conditions, i.e. m = 10 — | because 
the form of the bands involved in the transport of major- 
ity carriers (determining their transmission coefficient) is 
very robust under a slight change of the energy offset, see 
Figs. [3] and HJ_ 

From Fig. we clearly observe non-Ohmic behavior 
which is due to the rather complicated electronic struc- 
ture involved in the transmission probability. We ob- 
serve, for instance, that the absolute value of the current 
transmitted through a structure consisting of six layers 
CrAs is higher for all voltages than when transmitted 
through four layers. These TV characteristics indicate 
that spin-filtering should also be realizable with very thin 
structures of CrAs, see Fig. [l2j which might be easier to 
fabricate. For n=4, the minority current shows nonlin- 
earities which we attribute to resonant transport medi- 
ated by states which, in the bulk, give rise to the at bands 
discussed above. Nevertheless, it has to be kept in mind 
that the theoretical prediction is less reliable for very thin 
structures because (i) the modeling of the CrAs layers is 
based on the Hamiltonian of bulk ZB CrAs and (ii) ef- 
fects from the interface will become more important for 
thin layers. 



In Fig. [TT] we show the TV characteris- 
tics as obtained for different layers thicknesses, i.e. 
(GaAs) m /(CrAs)^/(GaAs) m , where I = 4,6,8,10 and 
m = 5. The results shown in Fig. 11 stem from simula- 



tions in which the number of layers of GaAs to the left 
and the right of CrAs was set to m — 5 and kept constant, 
i.e. the electric field across the CrAs layer at a given volt- 
age increases with decreasing layer thickness. This trend 
follows the actual physical trend within the device and 
that of a self-consistent model. Moreover, we note that 
the actual form of the computed TV characteristic is only 



V. DISCUSSION AND CONCLUSIONS 

We have performed a model study of transport in 
GaAs/CrAs heterostructures which is based on the as- 
sumption that sufficiently thin layers of fee CrAs can 
be grown lattice matched in between a GaAs substrate 
and capping layer. The bulk electronic structure of 
fee GaAs, fee CrAs, as well as lattice matched sin- 
gle [1,0,0] GaAs/CrAs heterointerfaces were calculated 
within LSDA and used to determine the band offset 
between the two materials. For a lattice constant of 
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Figure 11: (Color online) I-V characteristics for different layer 
thicknesses £ = 4, 6, 8, 10. 
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Figure 12: (Color online) Current spin polarization for differ- 
ent layer thicknesses £ = 4, 6, 8, 10. 

ai = 5.65 A we find that fee CrAs is a half-metal, 
with zero gap for one spin orientation (majority carri- 
ers) and a gap of 1.8 eV at the X point for the other 
(minority carriers). The computed band offset of about 
0.5 — 0.6 eV, however, aligns the gap region of minority 
CrAs with the central region of the top valence bands 
of GaAs. For spin-filtering therefore it is not important 



whether the sheet of fee CrAs lattice-matched to GaAs 
is half-metallic. 

The computed ab-initio spin-dependent electronic 
structures are downfolded onto TB models which are 
used to construct the effective Hamiltonian of the of n- 
GaAs/CrAs/n-GaAs heterostructures consisting of £ = 
4, 6, 8, 10 mono-layers of CrAs and to compute the cur- 
rent response within a non-equilibrium Green's function 
approach. Downfolding was constrained by the require- 
ment that As TB onsite parameters for a given spin ori- 
entation be constant throughout the system. We con- 
sider carrier injection from n-doped GaAs and our calcu- 
lations show efficient spin-filtering over a wide parameter 
range, in particular, regarding the precise band align- 
ment between the GaAs conduction band edge with the 
CrAs bands, temperature, and layer thickness. Spin- 
polarization of up to 99 percent is predicted within this 
model. 

A number of potential improvements to the present 
theoretical approach, such as a more realistic account of 
correlation, inclusion of the spin-orbit interaction, a self- 
consistent treatment of transport, etc., can readily be 
listed and be addressed in future studies. However the 
main problem currently lies in the fabrication of fee het- 
rerostructures containing layers of transition metal com- 
pounds, such as CrAs, MnAs, or VAs, and conventional 
fee semiconductors. For MnAs, apart from Mn delta- 
doped GaAs structures and strained fee MnAs quantum 
dots on GaAs, the growth in the fee apparently has not 
been achievable. Evidence for fee CrAs layers on GaAs 
substrates still seems to be controversial. It is hoped that 
these promising theoretical results regarding high spin- 
polarization encourage the materials growth and experi- 
mental physics community in the study of semiconductor 
heterostructures containing transition metal compounds. 
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Table I: tight-binding parameters for bulk CrAs majority and 
minority spin and the respective GaAs parameters. The anion 
onsite energies are indicated by the number 1 while the cation 
is labeled by the number 2. For further notations see^. 
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CrAs min 
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4.3837 


4.1719 


-11.8119 


1.0108 


5.5583 


-16.4426 




1.2291 


1.2291 
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